Novel features of nuclear forces and shell evolution in exotic nuclei 
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Novel simple properties of the monopole component of the effective nucleon-nucleon interaction are 
presented, leading to the so-called monopole-based universal interaction. Shell structures are shown 
to change as functions of N and Z consistently with experiments. Some key cases of this shell 
evolution are discussed, clarifying the effects of central and tensor forces. The validity of the present 
tensor force is examined in terms of the low-momentum interaction Wi ow k and the Qbox formalism. 
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Exotic nuclei provide us with new phenomena which 
are not found in stable nuclei. One of them is the evo- 
lution of shell structure as function of proton number 
(Z) or neutron number (N) [H, Q. The evolution ends 
up, in some cases, with the appearance of new magic 
numbers and/or the disappearance of conventional ones. 
As Z increases, there are more exotic isotopes between 
the /3-stability line and the drip line, creating a wider 
frontier. Most of such exotic nuclei are far inside the 
drip line, being well bound. The driving force behind 
the change in their structure should be the combination 
of the unbalanced Z/N ratio and the nuclear force [![. 
Thus, it is crucial to see robust basic features of the nu- 
clear force in exotic nuclei. We present in this Letter 
novel simple properties of the monopole component of 
shell-model interactions which can reproduce experimen- 
tal data. While the shell evolution due to the tensor force 
has been suggested in 0], we introduce here an interac- 
tion which includes the central force also, moving closer 
to the complete picture. 

We start with selected shell-model nucleon-nucleon 
(NN) interactions which are successful in describing ex- 
perimental data. These interactions were obtained based 
on so-called microscopic interactions, derived for example 
with the G-matrix approach [B|, 0] starting from a bare 
NN interaction and incorporating short-range repulsion 
and core polarization. In order to reproduce experimen- 
tal data, however, the microscopic interaction has to be 
modified empirically, as is the case for the families of the 
USD 0], KB3 § and GXPF1 [| interactions. We shall 
take p/-shell first, and analyze the GXPF1A interaction 
flpj ] and the corresponding G-matrix interaction Q. 

The monopole matrix element of a given two-body in- 
teraction, V, is defined as 

_ E k , k '(jkfk'\V\jkfk') 



where j denotes a single-particle orbit with k being its 
magnetic substate and ( • • \ V\ ■ ■ ) is the antisymmctrized 
two-body matrix element. The monopole component of 
V is written, for j^f, as J2j j> v m;j,j' n j n j'i where rij is 
the occupation number (operator) of orbit j 12]. The 
monopole component is nothing but the average over all 
orientations. It was introduced by Bansal and French 
[lH , while its relevance to the effective shell model inter- 
action was discussed by Poves and Zuker [8[. Recently, 
the monopole component of the spin-isospin interaction 
has been shown to modify even the magic structure in 
exotic nuclei [3(, and the specific and substantial role of 
the tensor force was shown in Ref. Note that u m ;i,j' 
is defined either with isospin, T=0 or 1, or in the proton- 
neutron scheme, while "j, j'" may be omitted for brevity. 

The importance of the monopole interaction for exotic 
nuclei originates in its linearity. As the orbit j' is occu- 
pied, the single-particle energy (SPE) of an orbit j, Cj, is 
changed by [12j . 

Aej = v m ,j. r n r . (2) 

For j'=gg/2 as an example, njr takes values up to 10. 
Thus, the effect of the monopole component can be mag- 
nified considerably. By moving along the nuclear chart, 
one can indeed change a particular ny substantially. This 
highlights the physics of exotic nuclei compared to that 
of stable nuclei, and it is of keen and urgent interest to 
clarify general and robust features of the monopole inter- 
action. At shell closures, the monopole component pro- 
duces effects according to Eq. © , whereas effects of other 
multipolc components vanish. The monopole component 
governs (spherical) SPEs on top of closed (sub-)shells. 
In open shell systems, its effects can be viewed through 
Eq. (J2]) as effective SPEs. As the surface deformation 
with low excitation energies is a Jahn- Teller effect, the 
SPEs are crucial for collectivity too. 
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FIG. 1: (color online) Monopole matrix elements of various forces for (a-d) pf and (e-h) sd shells. In (b,d,f,h), the tensor-force 
effect is subtracted from the others, and results from a Gaussian central force are shown. 



Figure Q] (a) shows w m; j t j> for isospin T=0 from the 
GXPF1A interaction, the G-matrix interaction Q and 
the tensor force in the p/-shell. The tensor force refers 
hereafter to the ir+p meson exchange force used in 
The orbits (j,f) are grouped as (/, /), (p, p) and (/, p). 
In Fig. [T](a), we find two distinct kinks in the tensor-force 
values for the (/, /) and the (p, p) groups, and the same 
kinks appear also in the GXPF1A and the G-matrix re- 
sults. Note that each kink is a consequence of the general 
rule suggested in The similarities are remarkable. To 
shed more light on this, in Fig. Q] (b) we subtract the 
tensor-force contribution from the GXPF1A and the G- 
matrix values. This results in almost flat curves. The (/, 
/) and (p, p) cases show almost the same values, while 
the (/, p) shows higher but still nearly flat values. This 



can be understood in terms of radial integral of the cen- 
tral force: in the former case the radial wave functions 
are the same between j and j' , while they are different in 
the latter. The flatness suggests a longer-range central 
force. In order to incorporate these features, we intro- 
duce a central Gaussian interaction as 

Vc = J2 /s,TP S ,Texp(~(r/ M ) 2 ), (3) 

S,T 

where S(T) means spin (isospin), P denotes the projec- 
tion operator onto the channels (S,T) with strength /, 
and r and fi are the internucleon distance and Gaussian 
parameter, respectively. Figure [1] (b) shows results ob- 
tained by /o,o=/i, 0=166 MeV and /i=1.0 fm. The agree- 
ment with GXPF1A is remarkable, considering the sim- 
plicity of the model. Thus, we can describe the monopole 
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FIG. 2: (color online) Diagrams for the Vmu interaction. 



component by two simple terms: the tensor force gener- 
ates "local" variations, while the Gaussian central force 
produces a flat "global" contribution. It is worth men- 
tioning that /x=1.0 fm is reasonable from the viewpoint 
of NN interaction, and deviations from it, including the 
zero-range limit, worsen the agreement. 

Figure [T] (c) shows D m 's for T=l. They are grouped 
for pairs of j=f and the rest. The former corresponds 
to the standard BCS-type pairing cases. We first stress 
that the basic scale is quite different between T=0 and 1: 
v m 's of GXPF1A are in the range -2.5 - -1 McV for T=0, 
whereas for T = 1 they are in the range -0.3 ~ 0.2 MeV. 
The sharp rise for j=j'=Pi/2 occurs in all three interac- 
tions as a characteristic fingerprint of the tensor force. 
Note that v m 's for the GXPF1A interaction (G-matrix) 
are mostly repulsive (attractive) for j^j' . A similar re- 
pulsive correction to G-matrix values occurs in the sd- 
shell as well. We subtract the tensor contribution as was 
done in Fig.[T](b), and show the result in Fig.[T](d) as well 
as those of the Gaussian central force with /o.i=0.6/o,o 
and /i,i=-0.8/o,o- The basic feature can be reproduced, 
apart from some deviations in the (/,/) cases, which may 
indicate stronger pairing correlations. 

Figure [1] (e-h) exhibits v m 's in the sd-shell, similar to 
what is shown in Figs. CD (a-d). The SDPF-M interac- 
tion [l3| is taken as the realistic interaction. All fea- 
tures discussed for the p/-shell are seen, and the tensor- 
subtracted values are reproduced by the same Gaussian 
central force even better for T=l. One sees repulsive 
corrections to w m 's from the G-matrix for T=l and j^=j' , 
similar to our findings in the pf shell. This correction 
is linked with the oxygen drip line, its origin has been a 
puzzle, but has recently been resolved 14 1 . 

Based on the above results, we introduce the monopole- 
based universal interaction, Vmu- As shown in Fig. [51 
Vmu consists of two terms. The first term is the Gaus- 
sian central force discussed so far, and should contain 
many complicated processes including multiple meson ex- 
changes. The second one is the tensor force comprised 
of 7r and p meson exchanges [2]. The Vmu interaction 
resembles Weinberg's original model for Chiral Pertur- 
bation theory [l^l, if one replaces Fig. [2] (a) by contact 
terms and (b) by the one-7r exchange potential. 

Figure [3] shows applications of Vmu to the shell evolu- 
tion assuming a filling configuration. Figure[3] (a) depicts 



neutron SPEs around A"=20 for Z=8~20. Starting from 
SDPF-M SPEs at Z=8, one sees the evolution of the 
7V=20 gap, in a basically consistent manner with other 
shell- model studies [3, EH • While the change is mono- 
tonic without the tensor force, the tensor force produces 
a sharp widening from Z=8 to 14, and then stabilizes the 
gap towards Z=20. It is worth mentioning that the nor- 
mal SPEs arise at Z=20, whereas at Z=8 the inversion 
between f 7 / 2 and p 3 / 2 occurs and d 3 / 2 is rather close to 
P3/2, leaving the major gap at A^=16. The central force 
lowers the neutron d 3 / 2 SPE more than the fj/ 2 SPE as 
protons occupy the sd-shell due to larger overlaps, yield- 
ing a wide 7V=20 gap at 40 Ca. The N=20 gap at Z~14 
is, however, largely due to the tensor force, and becomes 
smaller if protons are excited to d 3 / 2 . 

Figure [3] (b) shows proton SPEs for 68_78 Ni starting 
from empirical values 17f at N=40. The SPE of pi/ 2 is 
not known empirically, and is placed above p 3 / 2 by the 
energy difference predicted by the GXPF1A interaction. 
The orbit f 5 / 2 crosses p 3 / 2 at N=45 consistently with 
a recent experiment [18J, and the f 7/2-/5/2 splitting is 
reduced by 2 MeV from iV=40 to 50. For both, the tensor 
force plays crucial roles. Apart from certain differences, 
the trend is seen in other shell model results, e.g., [l9j |. 

Figure[3] (c) shows neutron SPEs relative to d 5 / 2 on top 
of 90 Zr- 100 Sn, starting from empirical values at Z=40 ob- 
tained by averaging with spectroscopic factors [20| ■ The 
lowering of 177/2 is remarkable 21] . If there were no 
tensor-force effects, 37/2 and h\\/ 2 do not repel, ending 
up with quite a different shell structure for 100 Sn, making 
this nucleus much softer. The closer spacing of g 7 / 2 and 



d 5/2 in 101 



Sn seems to be seen experimentally [221 ] . 
We now discuss whether the simple tensor force in Vmu 
can be explained microscopically or not. We take the 
AV8' interaction [23| and examine how the tensor force 
obtained by the spin-tensor decomposition changes in the 
following processes. We derive a low-momentum interac- 
tion Viowk [2J] and calculate shown in Fig. 0] 
(a,b) for T—0, 1, varying the cut-off parameter, A. For 
the usual value A=2.1 (fm" 1 ), the result is very close to 
the bare AV8' tensor force contribution. We then perform 
the 3rd-order Q-box calculation with folded diagram cor- 
rections Q, in order to include medium effects like core 
polarization. The result still resembles K m 's of the bare 
tensor part. Thus we can confirm that the treatments 
of the short-range correlation and the medium effects do 
not change much i> m 's of the tensor force. This near- 
independence may be interpreted in terms of specific and 
complicated angular momentum coupling in the tensor 
force. For instance, the 2nd-order perturbation by two 
tensor forces yield mainly a central force. For unusual 
values like A=l (fm -1 ), deviations arise, as expected. 

The central force depends strongly on A. For A=2.1 
(fm -1 ), z) m 's of the central part of Vi ow k are scattered 
around the values of Fig. [1] (b). This result is promising, 
but more studies are needed. 
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FIG. 4: (color online) Tensor forces in AV8' interacton, in 
low-momentum interactions obtained from AV8', and in the 
3rd-order Qbox interaction for (a) T=0 and (b) T=l. 



In summary, we have presented novel general prop- 
erties of the monopole interactions, and introduced the 
Vmu interaction consisting of simple central and tensor 
forces. The persistency of the bare tensor force is ex- 
amined by the Yi ow k and Qbox formalisms. The Vmu 
produces a variety of the shell evolution, connecting sta- 
ble and exotic nuclei, for instance, exotic Ne-Mg with 
40 Ca, 68 Ni with exotic 78 Ni, and 90 Zr with exotic 100 Sn. 
The shell structure appears to vary considerably in exotic 
nuclei. As Vmu has been introduced based on monopole 
properties, tests of its validity by full shell-model calcula- 
tions and possible refinements including multipolc com- 
ponents are of great interest, as well as more predictions 
by Vmu- 



We thank Dr. H. Grawe for valuable discussions. This 
work was supported in part by a Grant-in- Aid for Spe- 
cially Promoted Research (13002001). This work has 
been supported by the JSPS Core-to-Core Prgram EFES, 
and is a part of the RIKEN-CNS joint research project 
on large-scale nuclear-structure calculations. 



[1 

[2 

[3 

[5 

[6] 

[7] 

[8] 
[9] 

[10 

[11 
[12 

[13 
[14 
[15 
[16 

[17 
[18 
[19. 
[20 
[21 
[22 

[23 
[24 



Glasmacher, Prog. Part. Nucl. Phys. 60, 
-G. Porquet, Prog. Part. Nucl. Phys. 61, 



Rev. Lett. 87. 
Rev. Lett. 95, 
Brown, Nucl. 



082502 (2001). 
232502 (2005). 
Phys. A114, 241 



Kuo and E. Osnes, Phys. 



A. Gade and T. 
161 (2008). 
O. Sorlin and M 
602 (2008). 

T. Otsuka et al., Phys 
T. Otsuka et al., Phys 
T.T.S. Kuo and G.E. 
(1968). 

M. Hjorth- Jensen, T.T.S 
Repts. 261, 125 (1995). 

B. A. Brown and B.H. Wildenthal, Ann. Rev. Nucl. Part. 
Sci. 38, 29 (1988). 

A. Poves and A. Zuker, Phys. Rep. 70, 235 (1981). 

M. Honma et al. Phys. Rev. C 65, 061301(R) (2002); 

Phys. Rev. C 69, 034335 (2004). 

M. Honma et al, Eur. Phys. J. A 25, sOl, 499 (2005). 
R.K. Bansal and J.B. French, Phys. Lett. 11, 145 (1964). 
For j—j' for identical nucleons, slightly different forms 
are used. 

Y. Utsuno et al., Phys. Rev. C 60, 054315 (1999). 
T. Otsuka et al., Phys. Rev. Lett., submitted. 
S. Weinberg, Phys. Lett. B 251, 288 (1990). 
F. Nowacki and A. Poves, Phys. Rev. C 79, 014310 
(2009). 

H. Grawe et al, Eur. 
K.T. Flanagan et al., 
N.A. Smirnova et al, 



http://www.nndc.bnl.gov/ensdf/ 



Phys. J. A 25, sOl, 357 (2005). 
Phys. Rev. Lett., submitted. 
Phys. Rev. C 69, 044306 (2004). 



P. Federman and S. 
D. Seweryniak et al. 
R. Gryzywacz et al. 
B.S. Pudliner et al, 
S.K. Bogner, T.T.S. 
386, 1 (2003). 



Pittel, Phys. Lett. 69B, 385 (1977). 
Phys. Rev. Lett. 99, 022504 (2007); 
priv. com. 

Phys. Rev. C56, 1720 (1997). 
Kuo, and A. Schwenk, Phys. Rep. 



